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Abstract 

A natural consequence of the composite operator propagator-vertex description of deep 
inelastic scattering developed by the authors is that the anomalous suppression observed 
in the flavour singlet contribution to the first moment of the polarised proton structure 
function gf (the 'proton spin' problem) is not a special property of the proton structure 
but is a target independent effect which can be related to an anomalous suppression in 
the QCD topological susceptibility. In this paper, it is shown how this target independent 
mechanism can be tested in semi-inclusive deep inelastic scattering in which a pion or D 
meson carrying a large target energy fraction z is detected in the target fragmentation 
region. 
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1. Introduction 

The anomalous suppression of the first moment, F^, of the polarised proton structure 
function has been the focus of intense theoretical and experimental activity for nearly a 
decade. While it is now generally accepted that the key to understanding this effect is the 
existence of the chiral U{1) anomaly in the flavour singlet pseudovector channel, there are 
several detailed explanations reflecting different theoretical approaches to the description of 
deep inelastic scattering (DIS) and proton structure. In this paper, we review one of these 
- the composite operator propagator-vertex (CPV) description of deep inelastic scattering 
developed by us in a series of papers [1-4] - and show how one of its key predictions, the 
target independence of the suppression mechanism, can be tested in future semi-inclusive 
DIS experiments. 

The essence of our approach is the decomposition of structure function moments into 
the product of perturbative Wilson coefficients, non-perturbative but target-independent 
composite operator propagators, and vertex functions describing the coupling of these 
operators to the target nucleon. 

The vertices, which are defined to be 'IPF with respect to a chosen set of operators, 
encode all the information on the structure and properties of the target. They are non- 
perturbative and not directly calculable, and play the same role in our formalism as the 
parton densities in the conventional QCD parton model description of DIS. However, just 
as the parton densities have a universal character, being equally applicable to DIS or 
hadron-hadron scattering, these IPI vertices also have a more universal role, being related 
in favourable cases to low energy nucleon couplings such as g^NN etc. They provide an 
alternative, complementary, description of the nucleon state. 

However, the most important feature of the CPV formalism as far as the 'proton spin' 
problem is concerned is the separation of the composite operator propagator from the 
target-dependent vertex. This allows us to distinguish between generic non-perturbative 
properties of QCD and effects which are characteristic of the particular target. Our pro- 
posal is that the anomalous suppression in the flavour singlet contribution to the first 
moment of gi is of the first kind, viz. a generic, target-independent feature of QCD, re- 
lated to the chiral U{1) anomaly but not special to any particular hadron. In fact [1-3], we 
are able to relate the relevant propagator to a fundamental correlation function in QCD, 
viz. the topological susceptibility x(0)) show that the suppression in P^ is due to an 
anomalously small value of its first moment x'(0). To confirm this interpretation, we have 
evaluated x'(0) using QCD spectral sum rules[4], and have found a suppression in good 
quantitative agreement with the current data[5,6] on gf. 

The natural next step is to see whether this target-independent suppression mecha- 
nism can be tested directly, by studying the structure functions of other targets besides 
the proton and neutron. (Unfortunately, for these two targets, isospin invariance already 
implies target independence.) The obvious choice for an alternative target is the photon, 
whose structure function gj may be measured in two-photon processes at a sufficiently high- 
luminosity e+e" collider. However, this turns out to be an exceptional case, since there is 
a direct axial current - two photon coupling via the electromagnetic chiral U{1) anomaly. 
The first moment sum rule for g'J{k'^), as a function of the target photon virtuality A;^, 
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has been presented in refs.[7,8], together with estimates of the relevant cross-section asym- 
metries in polarised colliders. The dependence of gj{k'^) on the virtuality displays many 
interesting features: QiiO) is zero by electromagnetic current conservation [9]; its asymp- 
totic value for /c^ greater than the hadronic scale is essentially given by the electromagnetic 
anomaly coefficient, with logarithmic corrections governed by the gluonic anomaly; and its 
detailed dependence on k'^ as the various quark thresholds are crossed depends critically 
on the realisation of chiral symmetry in QCD. 

Direct DIS experiments on other hadronic targets are of course not feasible. We can 
nevertheless test our ideas in semi-inclusive DIS in an appropriate kinematic region where 
the reaction is well described in terms of deep inelastic photon scattering off a Reggeon 
(or more complicated exchanged object) with well-defined hadronic characteristics. In 
particular, using the target-independent suppression hypothesis, we are able to formulate 
predictions for ratios of cross section moments (related to moments of the Reggeon struc- 
ture functions) which are significantly and characteristically different from expectations 
based, like the EUis-Jaffe sum rule for g^, on the simple valence quark model or OZI rule. 
In particular, our target-independent mechanism should be clearly testable by comparing 
the ratios of cross section moments for the semi-inclusive reactions ep — > e7r~(D~)X and 
en en + (D0)X, in which a pion or D meson carrying a large target energy fraction z is 
detected in the target fragmentation region. 

The paper is organised as follows. In section 2, we review the most important features 
of the CPV method and its application to the polarised proton structure function, explain 
why it leads to a target-independent suppression, and compare our prediction for the first 
moment with the most recent SMC data. Then, in section 3, we show how by assuming 
target-independence and exploiting flavour SU (3) symmetry we can derive predictions for 
ratios of structure function moments for a variety of hadrons, including some which differ 
dramatically from results using the OZI rule. 

Semi-inclusive DIS is introduced in section 4, where we use a combination of symme- 
try and dynamical arguments to show that our predictions for ratios of structure function 
moments can be realised as ratios of cross section moments in a certain kinematical re- 
gion. In this region, the cross sections may be written in terms of Reggeon structure 
functions, where the exchanged Regge trajectory has the required SU (3) properties. We 
then compare these results with the more precise description of semi-inclusive DIS in terms 
of fracture functions [10], and relate the Reggeon structure function to the recently intro- 
duced extended fracture functions [11]. We conclude with a summary of our predictions for 
the most interesting ratios of semi-inclusive cross sections. 
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2. Target Independence and Composite Operator Propagator- Vertex Method 



The starting point is the sum rule for the first moment of the polarised structure 
function g^, viz. 

Jo 

= ^C^'ic^siQ')) (a^ + la') + lc!{a,{Q'))a^>{Q') (2.1) 

Here, a^, and a°((5^) are the form factors in the forward proton matrix elements of the 
renormalised axial current, i.e. 

{p,s\Al\p,s) = s^^a^ {p,s\Al\p,s) = s^^^a' (p, s|^°|p, s) = s^a°(Q2) 

(2.2) 

where and are the momentum and polarisation vector of the proton. The de- 
pendence of the singlet form factor follows from the renormalisation of the singlet current 
described below. The perturbative series Cf^ (q;s(<5^)) and Cf (^as{Q^)) are OPE coeffi- 
cients and are now both known to O(a^) [12-14]. 

Because of the chiral U{1) anomaly, the singlet current A'l is renormalised and mixes 
with the topological density. Defining the bare operators = ^q^ij.'jdq and Qb = 
a^^fj-i^PT^y-Q^^Q^^^ we have (for n/ flavours) 

Q = QB-^^{l-Z)d>^A% (2.3) 

where Z is a divergent renormalisation constant. The associated anomalous dimension 7 
was flrst calculated in ref.[15] and is now known to 3 loops[16]. Matrix elements of A'^ 
therefore have a non-trivial scale dependence governed by 7. In particular, 

(2.4) 

where t = In ^ . 

The anomalous Ward identities for composite operator propagators are 

d^{0\T A^ O\0) - 2nf{0\TQ O\0) = {6aO) (2.5) 

where O denotes an arbitrary composite operator and 5aO is its chiral variation. Notice 
that with these definitions of the renormalised composites, the combination ^5^^^° — 2n/Q^ 
appearing in the anomalous Ward identities is the same for the bare or renormalised 



3 



operators [17]. The possibility of making such a definition is a consequence of the Adler- 
Bardeen theorem. 

The sum rule (2.1) is derived using the OPE for two electromagnetic currents. The 
dominant contributions arise from the operators of lowest twist and, within this set, those 
of spin n contribute to the nth moment of the relevant structure function. Eq.(2.1) is the 
special case for odd parity operators of twist 2 and spin 1, viz. 



J^(?)J-(-g) ~ 



(2.6) 



It is at this point that our CPV method and the conventional parton model analysis 
of DIS diverge. In the parton model, the form factors are related to quark and gluon 
densities as follows[18]: 



= Au - Ad 
= Au + Ad 



2As 



a°(Q2)^AE-n/^4?^A^(g2) 



2tt 



(2.7) 



where AE = Au + Ad + As and Ag = jo (^?(^) + ^?(^))- There is a scheme 
ambiguity in these identifications, which relate four parton densities to just three measured 
quantities. The above definitions, in which AE is chosen to be scale invariant to all orders 
(this is possible because of the Adler-Bardeen theorem), are made in the AB factorisation 
scheme[18-20]. 




Fig. 1: The description of DIS in the parton model. The upper hatched blob denotes the perturbative 
QCD corrections related to the Wilson coefficients in the OPE. 

This standard approach is illustrated in Fig. 1. The upper hatched blob represents 
the perturbative QCD corrections contributing to the coefficient functions Cf^ , Cf in 
the OPE. The factorisation theorems show that these diagrams, with two quark (gluon) 
propagators, give the leading contribution to the amplitude for large , thus allowing the 
simple parton interpretation of q{x) and g{x,i) as the probability distributions for finding 
a quark or gluon with momentum fraction x in the target proton. 

The EUis-Jaffe sum rule makes the assumption that As and A^r are zero in the proton. 
This is equivalent to the OZI (Zweig) rule prediction a° = a®. The crux of the 'proton 
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spin' problem is to understand the origin of the OZI breaking revealed by the measurement 
of r^, which shows that is strongly suppressed relative to its OZI value. Our favoured 
explanation in the context of the parton model is that the OZI breaking is due overwhelm- 
ingly to the gluon density Ag in eq.(2.7). We expect the OZI rule to apply to the scale 
invariant quark densities, so that (in the AB scheme) AS = a^, while the scale dependent 
Ag{Q'^) compensates to produce an anomalously suppressed a'^(Q^). This would accord 
with the central conjecture of our rather different approach, described below, and has the 
virtue of providing a scale invariant meaning to the OZI rule in the presence of the chiral 
U{1) anomaly. 

In our approach, we again start from the OPE but instead factorise the resulting ma- 
trix elements into the product of composite operator propagators and vertex functions. 




Fig. 2: The description of DIS in the composite operator propagator-vertex method. The double line 
denotes the composite operator propagator and the lower cross-hatched blob the 'IPF vertex function. 

This is illustrated in Fig. 2. To do this, we first select a set of composite operators Oi 
appropriate (see below) to the physical situation and define vertices as 'IPI' with 

respect to this set. Technically, this is achieved by introducing sources for these operators 
in the QCD generating functional, then performing a Legendre transform to obtain an 
effective action r[Oj]. The IPI vertices are the functional derivatives of r[Oj]. A generic 
structure function sum rule then takes the form 

\x F{x,Q') = j;cj")(Q2)(0|T of OdO)ra,,, (2.8) 

i j 

where O^^' are the lowest twist, spin n, operators in the appropriate OPE with Cf the 
corresponding Wilson coefficients. 

This decomposition splits the structure function into three pieces - first, the Wilson 
coefficients cf^Q"^) which control the dependence and can be calculated in pcrtur- 
bative QCD; second, non-perturbative but target-independent QCD correlation functions 
(composite operator propag ators) (0|T of di\0); and third, a non-perturbative, target- 
dependent vertex functions ^Q.pp describing the coupling of the target proton to the com- 
posite operators of interest. The vertex functions cannot be calculated directly from first 

The presentation here is a slight over-simplification. In general, there is a distinction between the 

cases where the OPE operators Oj are included in the set Oi and where they are not. See refs.[2,3] for a 
complete account of the Zumino (partial Legendre) transform formalism and its application to the 'proton 
spin' problem. 
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principles. They encode the information on the nature of the proton state and play an 
analogous role to the parton distributions in the more conventional parton picture. 

One of the main advantages of our method is that some non-perturbative information 
which is generic to QCD, i.e. independent of the target, is factored off into the composite 
operator propagator. This allows us to distinguish between non-perturbative mechanisms 
which are generic to all QCD processes and those which are specific to a particular target. 
Our contention is that the anomalous suppression in the first moment of is of the first, 
target-independent, type. 

As emphasised in refs. [3,4,21,22], it is important to recognise that this decomposition 
of the matrix elements into products of propagators and proper vertices is exact, indepen- 
dent of the choice of the set of operators Oi. In particular, it is not necessary for Oi to be 
in any sense a complete set. All that happens if a different choice is made is that the ver- 
tices ^Q.pp themselves change, becoming 'IPI' with respect to a different set of composite 

fields. Of course, while any set of Oi may be chosen, some will be more convenient than 
others. Clearly, the set of operators should be as small as possible while still capturing the 
essential physics (i.e. they should encompass the relevant degrees of freedom) and indeed a 
good choice can result in vertices ^Oipp which are both RG invariant and closely related to 
low energy physical couplings, such as Qt^nn or gTj.^^[?,,2?,]. In this case, eq.(2.8) provides 
a rigorous relation between high DIS and low-energy meson-nucleon scattering. 

For the first moment sum rule for it is most convenient to use the chiral anomaly 
immediately to re-express a°((5^) in terms of the forward matrix element of the topological 
density Q, i.e. 

= ^2n/(p|g|p) (2.9) 

where the matrix element, which scales with the anomalous dimension 7, is evaluated at 
the scale (52.(2) 

Our set of operators Oi is then chosen to be the renormalised fiavour singlet pseu- 
doscalars Q and $5 where, up to a vital normalisation factor, the corresponding bare 
operator is $55 = The normalisation factor[3,4] is chosen such that in the ab- 

sence of the anomaly $5 would have the correct normalisation to couple with unit decay 
constant to the U{1) Goldstone boson which would exist in this limit. This is important 

This quantity may be evaluated directly in lattice QCD. See ref.[24] for a brief review of the current 
status of lattice evaluations. Note that in order to incorporate fully the effects of the anomaly, it is 
necessary [25] to use dynamical fermions. 

^^"^ To be precise, what is referred to here is the 'OZI limit' of QCD, defined in ref. [25] as the truncation of 
full QCD in which non-planar and quark-loop diagrams are retained, but diagrams in which the external 
currents are attached to distinct quark loops (so that there are purely gluonic intermediate states) are 
omitted. This is a more accurate approximation to full QCD than either the leading large 1/Nc limit, the 
quenched approximation (small Uf at fixed Nc) or the leading topological expansion (Nc ^ 00 at fixed 
rif/Nc. In the OZI limit, the U{1) anomaly is absent, as is meson-glueball mixing[26], and there is an 
extra U{1) Goldstone boson. Notice, however, that no approximation is used in deriving eq.(2.10). The 
OZI limit is used here purely as the motivation for choosing a particularly convenient normalisation for 
^5. 
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later in justifying the use of the OZI approximation for the vertex, which is then RG 
invariant. 

We then have 



^1 singlet = ^^2n/ Cf (^.(Q^)) 



{0\T Q Q\0)Tqpp + (0|TQ$5|0)r$ 



5PP 



(2.10) 



where the propagators are at zero momentum and the vertices (which in this equation have 
the external proton wave functions amputated) are IPI wrt Q and ^>5 only. 

The composite operator propagator in the first term is the zero-momentum limit of 
the QCD topological susceptibility x(^^)) viz. 

x(/c^) = J dxe'^-''i{Q\T Q{x) Q{0)\0) (2.11) 

The anomalous chiral Ward identities show that x(0) vanishes for QCD with massless 
quarks, in contrast to pure Yang-Mills theory where x(0) is non-zero. Furthermore, it can 
be shown[3,4] that the propagator (0|T Q $5|0) at zero momentum is simply the square 
root of the first moment of the topological susceptibility. We therefore find: 

n Singlet = ^ ^2n/ (Q^)) VxW T^,,, (2.12) 

The quantity y^x'JO) is not RG invariant and scales with the anomalous dimension 7. On 
the other hand, the proper vertex has been chosen specifically so as to be RG invariant. 
The renormalisation group properties of this decomposition are crucial to our resolution 
of the 'proton spin' problem. 

Our proposal (which is fully motivated in refs.[3,23] and supported by a range of low- 
energy phenomenology in the U{1) channel, such as rj' — > 77 decay) is that we should 
expect the source of OZI violations to lie in the RG non-invariant, and therefore anomaly- 
sensitive, terms, i.e. in x'(0).("') Since the anomalous suppression in is assigned to the 
composite operator propagator rather than the proper vertex, the suppression is a target 
independent property of QCD related to the chiral anomaly, not a special property of the 
proton structure. This immediately raises the question whether it is possible to test the 
mechanism by effectively performing DIS experiments on other hadronic targets. 

Our quantitative prediction then follows by using the OZI approximation for the 
vertex F$gpp and a QCD spectral sum rule estimate of the first moment of the topological 
susceptibility. (5) We find, for n/ = 3, 

Vx'iO) = 23.2 ± 2.4 MeV (2.13) 

^ ^ ^ ^ Q2 = iOGeV^ ^ ' 



'^^^ Notice that we are using RG non-irivariarice, i.e. dependence on the anomalous dimension 7, merely 
as an indicator of which quantities are sensitive to the anomaly and therefore likely to show OZI violations. 
An alternative suggestion, in which the suppression in a°((5^) is due directly to non-perturbative effects 
in 7 at low scales, was made in ref. [27]. This would also predict a target-independent suppression. 

The validity of this calculation has been criticised by IofFe[28,29] (see also ref.[30,31]), who asserts 
that the spectral sum rule technique cannot be applied to the ?7(1) channel because of problems with 
the optimisation scale and dependence on the strange quark mass. In ref. [32], we extend our analysis to 
include light quark masses and explain in detail why these criticisms are not valid. 
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This is a suppression of approximately a factor 0.6 relative to the OZI value f^/V^- 
Our final result is then 



a 



"(Q^^lOGeV^) = 0.35 ±0.05 (2.14) 



from which we deduce 



= 0.143 ±0.005 (2.15) 

Q2=ioGeV2 



This is to be compared with the EUis-Jaffe (OZI) prediction of — 0.58 ± 0.02 and the 
SMC experimental data[5] : 



= 0.136 ±0.013 ±0.009 ±0.005 (2.16) 

Q2 = 10GeV2 



where the last error is theoretical, related to the evolution. This gives 

a°(g2 = lOGeV^) = 0.28 ±0.16 (2.17) 

There is, however, a remaining uncertainty over the data related to the small x region. 
The SMC experiment is limited to measuring the region 0.003 < x < 0.7, and only a small 
estimated contribution of 0.0042 ±0.0016 is included in eq.(2.16) for the contribution to 
from the unmeasured range < x < 0.003. (The high x extrapolation is uncontroversial.) 
Recent fits[19,33] to the same data using a different extrapolation to the small x region, 
incorporating evolution of the parton distributions, suggest a much smaller central 

value for a° with larger errors, viz. a° 

(g2 = lOGeV^) = 0.10 Interestingly, these 
fits also suggest that AE = 0.45 ± 0.09, not too far from the OZI value. 

Very recently, new preliminary proton data has become available from SMC [6]. This 
gives 

[ dx g^{x; = lOGeV^) = 0.146 ± 0.006 ± 0.009 ± 0.005 (2.18) 

io.003 

The result for the entire first moment depends on how the extrapolation to the unmeasured 
small X region is performed. Using a simple Regge fit, SMC find = 0.149 ± 0.012 from 
which a° = 0.41 ±0.11, while using a small x fit using perturbative QCD evolution of the 
parton distributions they find Ff = 0.135 ± 0.016 which implies a° = 0.27 ± 0.15 (all at 
g2 = lOGeV^). 

Clearly, much more analysis, both theoretical and experimental, of the small x be- 
haviour of the polarised structure functions is required and studying this region will be an 
important goal of future polarised collider experiments at HERA. Nevertheless, the broad 
agreement with our prediction (2.14) is very encouraging and strongly suggests that our 
interpretation and explicit calculation of the topological susceptibility are correct. 
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3. The gi Sum Rule for Other Targets 



In this section, we consider the imphcations of the target-independent suppression 
mechanism for the structure functions of other hadrons, leaving aside temporarily the 
question of how this may be realised experimentally. 

Our basic predictionC^) is that for any hadron, the singlet form factor in eq.(2.1) can 
be substituted by its OZI value multiplied by a universal (target-independent) suppression 
factor s(Q^) determined, up to radiative corrections, by the anomalous suppression of the 
first moment of the topological susceptibility a/x'(0). For example, for a hadron containing 
only u and d quarks, the OZI relation is simply a° — a^, so we would predict: 

r? = ^C^'iasm (a^ + ^(1 + 4.)a«) (3.1) 

where 

Since s is target independent, we can use the value measured for the proton to deduce Fi 
for any other hadron target simply from the flavour non-singlet form factors, which obey 
relations from flavour SU{?>) symmetry. From our spectral sum rule estimate of a/ x'(0), 
we flnd s ~ 0.66 at Q"^ = lOGeV^, while the central value of the SMC result (2.17) gives 
s ~ 0.55. (We use the experimental data taken directly from SMC, ref.[5] in this section.) 

The form factors for a hadron B are given by the matrix elements of the flavour octet 
axial currents. The SU (3) properties are summarised by 



{B\4'ly\B) = {i^ii-i^iim 



fB yB jB yB 




(p^|A(^)|p^) (3.3) 



Here, p indicates the SU (3) representation while /, I3 and Y are the isospin and hyper- 
charge quantum numbers. The term |p^) is a reduced matrix element, while the 
other factors are SU{2) and S'C/(3) Clebsch- Gordon coefficients [36]. 

If we now take the hadron B to be in the 10 representation, then since 

10 X 10 = 1 + 8 + 27 + 64 (3.4) 



(6) The analo gous prediction in the parton model would be to assume that (in the AB scheme) the RG 
invariant AS would take its OZI value, ie. AS ~ where AE^,a; is the sum of the valence quark 

densities. These can be distinguished from the OZI-violating sea quark densities in semi-inclusive DIS in 

the current fragmentation region[34,35]. The gluon contribution AG = nj ""2^ Ag(Q-^) would then be 
given by 

AG = (1-5(Q')) AS^„, 

where s(Q^) is simply (3.2) with the Wilson coefficients omitted. This has the correct scaling property 
and ensures a target independent suppression in Pi. 
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the matrix element of the (octet) currents contain just one reduced matrix element. This 
is in contrast to the case of B in the octet representation, as for the proton or neutron, 
which would involve an F/D ratio arising from the two reduced matrix elements in the 
decomposition 8x8 = 1 + 8 + 8 + 10 + 10 + 27. This is an important simplification, as it 
means that the ratio of Fi for decuplet states can be predicted as a simple group-theoretic 
number, up to the dynamical suppression factor s. 

For example, for the A"'""'", the matrix elements of the currents are 

(A++|y4j|A++) = y5"(10|yl(«)|10) {A++\Al\A++) = ^{10\A^^^\10) (3.5) 

evaluating the relevant Clebsch-Gordon coefficients. Similar results hold for the A~. Tak- 
ing the ratio to eliminate the common reduced matrix element (10|A*^^) |10), we find the 
following result for the ratio of the first moment of the polarised structure functions gi for 
the A++ and A-: 



rf^ ^ V^ + Vll)Vl(l + 4^) ^ 2s + 2 



The OZI (c.f. EUis-Jaffe) prediction is given by setting s = 1, i.e. Ff /Ff = 4.(7) How- 
ever, substituting a suppression factor of s ~ 0.66 gives a much larger ratio Ff^ /rf' ~ 
10, while the experimental factor s ~ 0.55 would give an even larger value, indicating a 
near complete suppression of Ff^ . 

We would therefore expect to find a quite spectacular deviation from the quark model 
expectation for this ratio of structure function moments. We can also show (footnote 

E++ E° 

(7)) that the same result is obtained for the ratio F^^ /^i" for the charmed baryons 
E++ = uuc and E[! = ddc. Of course, these examples have been specially selected (because 
of the 2s — 1 factor) to show a particularly striking difference from the simple valence quark 
model predictions. However, as we shall see in section 4, they are also the examples which 
can be dynamically isolated in semi- inclusive DIS. 

Although these are the most interesting, other ratios of structure function moments 
can be easily calculated by the same method. The most obvious is the proton-neutron 



Alternatively, this result can be simply obtained in the valence quark model as follows. Using the 
quark charges and neglecting radiative corrections, we have 

Ti = ^(4Au + Ad + As) 

and so 

rf = |au(a++) rf- = iAd(A-) 

3 6 



With the (isospin) assumption Au(A++) = Ad(A~), we immediately find T^^^/r^ = 4. 

S++ s° 

The corresponding result for the ratio Fj^ /Pj^ follows from quark counting in the same way, 
assuming Au(E^"'") = Ad(E°) and treating the heavy c quark as a spectator. 
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ratio which, as noted above, does not simply give a group theoretic number but depends 
also on the F/D ratio. In this case, the OZI prediction is 



~ 4 - e>F/D 

while including the anomalous suppression factor, we find 

r? _ 2s - 1 - 3(2s + 
ff ~ 2s + 2 - QsF/D 

This complements the Bjorken sum rule for the difference V\ —T^^ viz. 



(3.7) 



(3.8) 



rf-r- = ^-C^\a,{Q^)) gA (3.9) 
where qa — dp- The neutron structure function is [5] 

rr|g2=ioGeV^ = -0.046 ±0.021 (3.10) 
so that the experimental result for the ratio is 



J- 1 



= -2.96 ±1.39 (3.11) 

Q2 = 10GeV2 



This is to be compared with the OZI result F^/F" = —7.6 ± 1.4, where we have used 
F/D = 0.575 ± 0.016 in eq.(3.7), and with the prediction from our modified formula (3.8) 
which gives central values T^/T^ = —3.5 for s ~ 0.66 and Ff/F" = —2.9 for s ~ 0.55. 
Of course, this only confirms that the suppression in the singlet form factor aP{Q'^) is the 
same for the proton and neutron, as expected by isospin symmetry and confirmed by the 
experimental validity of the Bjorken sum rule. 

As a final example, we quote the corresponding results for hadrons containing the 
strange quark. For the octet E+ and S~, the OZI rule gives respectively aP = |(3a^ — a^) 
and a° = — |(3a'^ + a^), so that our prediction is 

= ^^^rK(Q')) ((1 + 2s)a' + 1(1- 2s)a') (3-12) 



while 

r 



r = ^C^'ic^siQ')) ((1 - 2s)a' ± 1(1 - 2s)a«) (3.13) 
A similar group theoretic calculation then gives 

Ff^ 2s - 1 - 3(2s ± l)F/i:» 



Ff " 2s - 1 - 3(2s - 1)F/D 
11 



(3.14) 



with the valence quark model prediction again being recovered by setting s = 1. 

On the other hand, for the 2*+ and E*~ in the decuplet, the F/D ratio is absent and 
instead we find the simple ratio 

2s + 1 

We shall see to what extent these predictions may be tested in the next section. 



4. Semi-Inclusive Deep Inelastic Scattering 

Of course, it is not possible to measure structure functions directly for baryonic targets 
such as the A or Sc. However, it is possible to test the ideas in the previous sections 
in semi-inclusive DIS reactions eN ehX, where h is a detected hadron in the target 
fragmentation region. 

In our case, we are interested in reactions where N is a nucleon target and the detected 
hadron h is, for the interesting cases described above, a pion or D meson. The electron 
can of course represent any lepton. There are distinct contributions to this process from 
the current and target fragmentation regions, which we require to be clearly distinguished 
kinematically. A large rapidity gap is therefore required between h and the inclusive 
hadrons X, with h in the target fragmentation region, h is also required to carry a large 
target energy fraction z (defined below). 

4.1 Single Reggeon Exchange Model 




h(P2) 

Fig. 3 : The semi-inclusive DIS reaction eN ehX in the target fragmentation region with z ^ 1 modelled 
by the exchange of a Reggeon B. 

In this kinematical regime, the process may be modelled as shown in Fig. 3, in which 
the exchanged object is a Reggeon B with well-defined SU{?>) quantum numbers. With a 
polarised beam and target, this kinematics allows us to measure the structure function gf 
of the exchanged Regge trajectory B. 



12 



This is analogous to the measurement of the photon structure function gl in polarised 
e"''e~ scattering in the DIS rcgion[7,8] (Fig. 4) or the pomeron structure function F2 in 
diffractive ep scattering[37] (Fig. 5). 

e 



e 



e 




e 



Fig. 4 : The deep inelastic two-photon process in polarised e+e scattering used to measure the structure 
function of the photon. 

e 



e 




Fig. 5 : The diffractive exchange process in ep scattering used to measure the structure function of 
the pomeron. 

Since the results of section 3 depend solely on the SU{2>) properties of the baryon 
they will still hold here despite the fact that B is interpreted as a Reggeon. Indeed, it is not 
even necessary (see section 4.2) to assume that the exchanged object is a single Reggeon 
- our final predictions for cross section ratios hold independently of the dynamical nature 
of the exchanged object, which could in principle be a multi-Regge exchange or more 
complicated structure, provided the SUiZ) properties are correct. 

If the target is a nucleon N and the detected hadron is an octet meson (tt), SU{Z) 
symmetry shows that B belongs to a representation on the rhs of 

8x8=1 + 8 + 8 + 10 + 10 + 27 (4.1) 

Since the 27 requires a 5-quark state, it is a good dynamical approximation at sufficiently 
large z that the 10 dominates the 27. However, there is no such argument for 8 domi- 
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nance over the 10. To isolate a unique representation for B, we must therefore choose a 
combination of N and h giving I^, Y quantum numbers for B which appear in the 10 but 
not in the 8. This is satisfied by the A^"*" and A~, as in section 3. The required ratio 
of first moments Vf- /T^ — ff^f , where now A"'""'" and A~ are Reggeons, is therefore 
obtained by comparing the reactions ep — > e7r~X and en — > ctt+X. 

These symmetry considerations are easily pictured by drawing valence quark diagrams 
for the WaB vertex. Fig. 6 shows the quark structure of the ep — > c7r~X reaction, while 
the corresponding 5-quark, 27 represenation, exchange is shown in Fig. 7. 




(10) 



Fig. 6 : Quark diagram for the NhB vertex in the reaction ep — > evr X with the Reggeon B in the 10 
representation. 




A++ (27) 



Fig. 7 : Quark diagram for the NhB vertex in the reaction ep 
representation. 



eTT X with the Reggeon B in the 27 



As in section 3, we find the same ratio holds for the moments " /F]^ °, which can be 
realised (Fig. 8) in reactions in which a D meson is detected by comparing ep — > eD~X and 
en — s> eD*^X. To justify the assumption made there of treating the c quark as a spectator, 
we must select events in which there is no charmed jet in the current fragmentation region. 




Fig. 8 : Quark diagram for the NhB vertex in the reaction ep eD X where the Reggeon B has quantum 
numbers of the S J"*" . 

Trajectories with the quantum numbers of the E"*" and S~ would be exchanged in the 
reactions ep — > eK°X and en — > eK+X (substituting an s quark for the c quark in Fig. 8). 
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However, with these reactions there are two possibihties for the exchanged trajectory, with 
either the E"*" (E~) in the 8 or the E*"*" (E*~) in the 10 being possible (in addition to the 
Zweig suppressed 27 contribution). As we saw in section 3, these give quite different ratios 
of structure function moments. Unfortunately, it would be difficult to distinguish these 
possibilities experimentally. In particular, there is no dynamical justification for assuming 
10 dominance over the 8. This is why we had to choose the A"'""'' or (E°) 

trajectories to obtain a clear test of the predictions of section 3. 

The dynamics of these semi-inclusive reactions in the large 2, target fragmentation 
region (Fig. 5) can be deduced by analogy with the photon structure function or diffractive 
pomeron exchange processes. The first moment of the polarised structure function for the 
Reggeon B is found from the polarisation asymmetry of the differential cross section in the 
target fragmentation region, i.e. 

f^~^ f] \ ^target v A.'wfy'^ 

Here, x = xb = z = ^ so that 1 - z = ^, y = t = -(p2 - p'o? = -/c^ 

and Yp = ^{2 — y). This kinematics is described further in Appendix A. 

If we now take the ratio of eq.(4.2) for the two reactions ep en~X. and en e7r"'"X 
(or ep — > eD~X and en — > eD^X), the factorised Reggeon emission factor Af{z, t) cancels 
out, leaving the ratio of structure function moments Fi predicted in section 3 to be given 
simply by the ratio of the cross section moments. Our final predictions for the cross section 
ratios are summarised in section 4.3. 

As well as predicting the ratios, which contain the essential physics we wish to test, 
we should also consider the absolute size of the relevant cross section asymmetries. In 
particular, we must check that the cross sections do not fall off too quickly as z approaches 
1 for our predictions to be seen clearly in the data. Returning to eq.(4.2), and making the 
ansatz that the Reggeon emission factor /S.f{z,t) appropriate to the polarised amplitude 
is the same as that for the unpolarised case, we expect 

A/(^,t) ~ F(t)(l-^)i-2«s(*) (4.3) 

where asit) is the Regge trajectory for the B. For the reactions of interest, the relevant 
trajectory is the A, for which a/^{t) ~ 0.0 + 0.9t, with t in GeV^. For the relevant 
experimental condition of t close to zero (corresponding to a small scattering angle 6 lab 
of h relative to the incident nucleon in a collider experiment), the cross section moment 
therefore falls off only as 1 — 2;. 

4.2 Fracture Functions 

Within the general framework of the parton model, the appropriate description of 
events in the target fragmentation region in semi-inclusive DIS is with fracture functions, 
introduced in ref.[10]. In this section, we show briefly how the fracture function description 
gives a more rigorous foundation to the results given in the previous section in terms of 
Reggeon structure functions. 
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Fig. 9 : The contribution to semi-inclusive DIS from the current fragmentation region. 




Fig. 10 : The contribution to semi-inclusive DIS from the target fragmentation region. 

The two distinct contributions to semi-inclusive DIS from the current and target 
fragmentation regions are shown in Figs. 9 and 10. The current fragmentation events 
are described by parton fragmentation functions D'l{z;Q'^), where i denotes the parton, 
while the target fragmentation events are described by fracture functions M/*^(x, 2; Q^) 
representing the joint probability distribution for producing a parton with momentum 
fraction x and a detected hadron h carrying energy fraction z from a nucleon N. 

The differential cross section for polarised, semi-inclusive DIS has been given in 
refs. [38,39] , including NLO corrections. (The equivalent results for the unpolarised case 
were calculated in ref.[40].) For our purposes here, we just quote the lowest order result: 



X 



dAa 
dxdydz 



Yp A-Ko^ 



1 



X 



-Aqi{x;Q^)D'^ 



1 



X 



(4.4) 

where we have expressed the result in terms of the variable z = (see Appendix A). Here, 

Aqi{x;Q'^) and AMj^^ {x, z; Q^) are the polarisation asymmetries of the parton densities 
and fracture functions respectively. Restricting to the target fragmentation region, we 
simply have: 



X 



dxdydz 2 s 



J2el AMr{x,z-Q') 



(4.5) 
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In the kinematical region 2; ~ 1 where the dominant process can be modelled as 
Rcggcon exchange (Fig. 5), we can compare this expression to eq.(4.2). In this limit, 
therefore, we can relate the Reggeon structure function to the fracture function as follows: 

e- ' dx AM/^^(x, z; Q^) =^ J dt Af{z, t) dxs gfixB, t; Q^) (4.6) 

This is just the first moment of the more general relation 

AM/^^(x,2;g2) = ! dt Af{z,t) gf(j^,t;Q') (4.7) 
^ — ' 2~1 / i — z \i — z / 

This relation expresses the Reggeon structure function gf in terms of its partonic 
constituents, as described by the fracture function AM^^^ . We therefore see that the 
fracture function measures the parton distribution of the exchanged object [10]. Indeed, this 
interpretation is more general than the particular relation (4.7), since the fracture function 
description is not dependent on a particular model (such as a single Rcgge trajectory) for 
the exchanged object. For example, if the process is modelled by multi-Regge exchange, 
the rhs of eq.(4.7) would comprise a sum over the Reggeons. 

We can take this identification a stage further by considering the extended fracture 
functions M/*^(a;, z,t; Q'^) introduced recently in ref.[ll]. These are defined such that 

Ml^'^{x,z;Q^) = / dtM^^{x,z,t;Q^) (4.8) 

where t = —{vi~v''i)^- Just as in the integrals of the Reggeon emission factors, the upper 
limit of the t integration is not precisely specified, with the physical results for large 
being independent of the exact choice to the required order. These extended fracture func- 
tions have a number of important features[ll], notably a much simpler, homogeneous, RG 
evolution equation. They also have an interesting interpretation in terms of spacelike cut 
vertices, whose RG properties are known to be determined by the anomalous dimensions 
of appropriate local operators. 

These extended fracture functions allow us to remove the t integration in eq.(4.7), so 
we can finally write, to leading order, 

J]efAM/^^(a;,^,t;Q^) = F(t)(l - ^)-2«- W (-^, t; Q^] (4.9) 

where we have substituted eq.(4.3) for Af{z^ t). The NLO corrections to the rhs of eq.(4.9) 
can be read off from refs. [38,39]. In fact, this should be taken as the definition of the 
'Reggeon structure function' in the single Regge exchange approximation. 

Since i3 is a Reggeon, we cannot express the structure function gf in terms of Wilson 
coefficients and operator matrix elements as for single particle structure functions such 
as gf. Nevertheless, the moments of gf should inherit a RG scaling dependendence on 
the anomalous dimension of the appropriate OPE operator, e.g. for the flavour singlet 
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first moment. For consistency, therefore, we would require the extended fracture function 
AM/*^(x, z,t; Q"^) to satisfy a homogeneous RG evolution equation. This is borne out by 
the results of ref.[ll], where it is shown that 



(4.10) 

where AP^- is the usual DGLAP evolution kernel. 

A satisfying picture therefore emerges, in which the results of section 4.1 are confirmed 
and placed in a broader theoretical framework for the description of semi-inclusive DIS. 



4.3 Predictions 



Of course, the ratio (3.6) is only obtained in the limit as z approaches 1, where the 
reaction eN — > ehX is dominated by the process in which most of the target energy is 
carried through into the final state h by a single quark (see Figs. 6-8). 

At the opposite extreme, for z approaching 0, the detected hadron carries only a 
small fraction of the target nucleon energy and has no special status compared to the 
other inclusive hadrons X. In this limit, therefore, the ratio of cross section moments (4.2) 
for ep — > e7r~X and en e7r"'"X is simply the ratio of the structure function moments for 
the proton and neutron, i.e. /FJ as given in section 3. The same result would hold for 
the ratio of ep — > eD~X and en — > eD°X. 

We therefore predict the following results for the ratios of the differential cross section 
moments ~' dx x f/aydZt ■ 



en^ e7r+(D0)X 
ep e7r-(D-)X 



2s -1 
2s + 2 
2s + 2 



6sF/D 



2s-l-3{2s + l)F/D 



(z^l) 
{z^O) 



(4.11) 



Between these limits, we can only interpolate. We therefore expect a plot of the ratios 
of ^ dx X '^^'^y^^^f over the range < 2 < 1 for en ^ e7r+(D°)X and ep e7r~(D~)X 
to look like the sketch in Fig. 11, where the solid line shows the ratios predicted by (4.11) 
with s ~ 0.66, contrasted with the ratios predicted by the OZI rule, i.e. s = 1 (dotted 
line) . 
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0.30 




-0.30 

Fig. 11 : Sketch of the cross section moment ratios for en e7r+(D'')X and ep eir^ (D^)X, interpolating 
between the Hmits z ^ and z ^ 1. The dotted line shows the OZI prediction and the solid line our 
prediction based on the target-independent suppression mechanism (with s ~ 0.66). 

The difference between the OZI (or valence quark model) expectations and these 
predictions based on our target-independent interpretation of the 'proton spin' data is 
therefore quite dramatic, and should give a clear experimental signal. 

Since our proposed experiment requires particle identification in the target fragmenta- 
tion region, it is difficult to do at a polarised fixed-target experiment such as COMPASS [41] 
at CERN, which is better suited to studying semi- inclusive processes in the current frag- 
mentation region. A better option is a polarised ep collider, such as HERA [42]. Testing 
our predictions requires comparision of proton and neutron data, which can be extracted 
from experiments with polarised deuterons replacing the protons in the collider. 
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Appendix A 



A number of different definitions of the variable 'z' in semi-inclusive reactions are 
used in the literature. Here, we describe the relation between our notation and that used 
elsewhere, and present a number of useful kinematical results. 

In refs. [38-40], the kinematics is described in the CM frame of the virtual photon and 
target nucleon. Let Eh and be the energy of the detected hadron and nucleon in this 
frame and 9 be the angle between the corresponding momenta. With the definition [40] 
f = 1(1 — cos 6*), we see that the target fragmentation region is characterised by f ~ 1, 
while the current fragmentation region is v ~ 0. The hadron energy fraction variable used 
by ref.[40] is then 

= E^l^ ^^-^^ 
In contrast, the corresponding variable used in refs. [34,35] is 

P2-q 

The relation is 

Notice [40] that these two variables are approximately equal in the current fragmentation 
region but differ substantially in the target fragmentation region, where is small. 

In terms of the variable t = — (p2 — ^2)^; which in the model of Fig. 5 is the invariant 
spacelike momentum —k^ of the exchanged Reggeon, we have 

xt 

SO that at fixed x, in the target fragmentation region, is simply a measure of t. The 
angle 9 (assuming t is small compared to Q^) is given by 

x{\-x)Q^ ^ ' 

Our preferred variable z = can be expressed in this frame as 



= (A6) 



= ^ 

^ En Q2 V- -^-(^) ■ -\Q 

so for relatively small t, z is simply given by the ratio of the detected hadron energy to the 
target nucleon energy in the photon-nucleon CM frame. The required kinematical region, 
where the semi-inclusive reaction is well approximated by the Reggeon exchange diagram 
and our prediction for the cross section moment ratios holds, is therefore v ~ 1 and z 
approaching 1. 
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